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Abstract 

We propose a novel stochastic radio resource allocation strategy that achieves long-term fairness considering 
backhaul and air-interface capacity limitations. The base station is considered to be only powered with a finite 
battery that is recharged by an energy harvesting source. Such energy harvesting is also taken into account in 
the proposed resource allocation strategy. This technical scenario can be found in remote rural areas where the 
backhaul connection is very limited and the base stations are fed with solar panels of reduced size. Our results 
show that the proposed scheme achieves higher fairness among the users and, in some cases, a higher sum-rate 
compared with the well-known proportional fair scheduler. 

I. Introduction 

We consider in this paper a downlink (DL) radio resource allocation strategy for a system with limited 
hackhaul capacity and where the base station (BS) is equipped with a finite battery recharged by an energy 
harvesting source. Although backhaul availability has been taken for granted in conventional systems, backhaul 
is, in general, a limited resource. This is the case, for instance, of the deployment planned in the European project 
TUCAN3G project (http://www.ict-tucan3g.eu) in which 3G femtocells empowered by solar panels of limited 
size and connected to the core network through a limited capacity WiFi-LD backhaul are used to provide 3G 
connectivity to users located in rural remote locations in Peru. Such limited WiFi-LD connections already exist 
currently and are used basically to provide remote health services. The final social objective of the project is to 
contribute to the economical development of such areas through the provision of communication services to the 
general users beyond the current limited health services. 

The backhaul capacity limitation can be introduced in the resource allocation problem by imposing a maximum 
instantaneous aggregate traffic rate constraint [1], [2], [3]. However, in real deployments, the backhaul capacity 
can only be measured in average terms. In addition, limiting the sum-rate instantaneously at each specific 
scheduling period may hamper the performance of the system in terms of the achievable long-term rates. In 
these circumstances, it seems less limiting to use high data rates in the access network whenever the channel 
conditions allow (even using greater instantaneous values than the average constraint imposed by the backhaul) 
provided that the average backhaul rate constraint is met when averaging the traffic served in several scheduling 
periods. Note that the backhaul constraint in terms of average traffic is suitable if we assume that queues are 
implemented at the entrance of the access network. 

In addition to the backhaul limitation, the energy available at the BS may be a very limited resource as well. 
If the BS is only powered with batteries (as may happen in rural environments, for instance in the TUCAN3G 
deployments described above), then the battery status as well as the harvesting capabilities (if any) should also 
be considered in the scheduling strategy explicitly if we want to optimize the performance subject to the energy 
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limitations. In the scenario that we will consider in this paper, the BS will he powered only hy a limited battery 
and an energy harvesting device, e.g., solar panels, which will allow recharging the battery [4]. 

In this paper we propose a long term fairness scheduler that considers a long-term backhaul constraint, the 
battery status of the BS, and the energy that it is being harvested. When there is no reason for treating flexible 
service rate users differently, the maximin criterion is a meaningful scheduling approach [5]. This approach 
maximizes, at each scheduling period, the minimum of the throughputs of the users. Essentially our goal is to 
provide a balanced long-term rate to a set of users. In addition, and as commented previously, the scheduler will 
consider explicitly the energy and backhaul constraints, while taking advantage in an opportunistic way of the 
instantaneous good wireless channel conditions. 

The remainder of this paper is organized as follows. In Section II we describe the system model. Section III 
presents the resource allocation strategy developed in the paper. The numerical evaluation, presented in Section 
IV, has been carried out using realistic models taken from remote rural locations in the forest in Peru. Finally, 
conclusions are drawn in Section V. 

II. System Model 

A. System Description 

Let us consider a DL scenario composed of a single BS and several users. We assume that the system is based 
on WCDMA technology [6] and two different types of users coexist: voice users and data users. Let us denote 
the set of voice and data users by ICy and ICd, respectively. We assume that voice users request a fixed service 
rate whereas data users request a flexible service rate. 

Users in WCDMA are multiplexed using codes [6]. We assume that the network operator has already reserved 
a set of codes for the voice users and the remaining codes are to be allocated among the data users. Thus, the 
amount of available codes in each set is known and fixed at the BS. 

The BS is solely powered with a battery and an energy harvesting source. The energy harvesting source allows 
the BS to collect energy from the environment and recharge the battery (for example, solar panels). This is 
especially important in rural areas, where the access to the power grid may be impossible or too expensive. We 
consider that only causal information is available for the resource allocation strategy, i.e., only information of the 
past and current harvesting collections and battery dynamics will be available to execute the scheduling strategy 
at each particular scheduling period, yielding to an online approach. 

One of the novelties of this work is that we account for a maximum backhaul rate constraint. However, instead 
of limiting the instantaneous access network data rates as the maximum flow allowed by the backhaul, as in 
[1], [2], and [3], we limit the average throughput served by the access network. That means that we allow the 
instantaneous rate in the access wireless links to surpass the backhaul limitation. This can be done whenever 
we have queues at the entrance of the access network and such queues are stable (which, in fact, is guaranteed 
by imposing that the average aggregated rate is not higher than the backhaul capacity). Moreover, as it usually 
happens with backhaul connections in real deployments, it is not possible in the transport layer to differentiate 
among different users with the same type of service, for example in DiffServ-based networks [7]. Accordingly, 
we consider in the resource allocation problem in the access network that the backhaul capacity is equally divided 
among the users with the same type of service. Fig. 1 presents the system architecture of the target rural scenario. 

B. Power Consumption Model and Battery Dynamics 

In this subsection, we introduce the power consumption and the battery model considered in this paper. The 
overall power consumption at the BS is modeled as the addition of the radiated power, which is divided into the 
power devoted to pilot channels (Tcpich assumed to be fixed) and the power consumed by the traffic channels 
and a fixed power consumed by the electronics of the BS {Pc), where t denotes the scheduling period. 
The model considered for the last term is based on [8] and includes the power consumption of the radio frequency 
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(RF) chains, the hasehand power consumption, and the consumption of the cooling systems. The maximum traffic 
power will he forced to depend on the current battery level of the BS, as it will he apparent later. 

Given that, the overall energy consumption hy the BS during the f-th scheduling period is 

E{t)^Ts-{PcPlCH + PBs{t) + Pc), yt, ( 1 ) 

where is the duration of the scheduling period. We consider that the amount of power that can he used for 
traffic services is limifed (due fo physical consfrainfs of fhe amplifiers of fhe BS equipmenf). Such maximum 
fraffic power will he denofed as P^^^, so Ppsit) < 

Lef B{t) he fhe energy stored af fhe haffery of fhe BS at the beginning of the scheduling period t. Then at 
period t + 1, the battery level is updated in general as 

B{t + l) = f{B{t),E{t),H{t)), yt, (2) 

where H{t) is the energy harvested in Joules during the scheduling period t and the function /(•) : M+ x M+ x 
M+ —>• M+ depends upon the battery dynamics, such as storage efficiency and memory effecfs. A common practice 
is to consider the following battery update: 

B{t + l) = {B{t)-E{t) + H{t))f^\ yt, (3) 

where (x)q is the projection of x onto the interval [a,b], i.e., (x)^ = min{max{a, x}, 6}, which accounts for 
possible battery overflows and assures fhaf fhe baftery levels are non-negalive, and is the battery capacity. 

Notice that the whole harvesting collected during period t is assumed to be available in the battery at the end of 
the period for simplicity. In general, the total energy consumed by the BS during one period will be limited by 
a function of the current battery level as 

Ts • (PcPICH + PBs{t) + Pc) < giPit)), yt, (4) 

where the function g{-) is defined as g{B{t)) = min{Ts (Pcpich + P^^^ + Pc), u>{B{t))}, and w{-) : M+ —)• M+ 
a generic confinuous increasing funcfion fhat fulfills fhaf w{B{t)) < B{t), Vf. For example, if all fhe haffery is 
allowed to be spent during one particular epoch, then w{B{t)) = B{t). Nevertheless, the approach followed in 
this paper is to limit the amount the battery that can be used in one particular period in order to use the energy 
in a more conservative way. Thus, we consider that only a given fraction of the battery is allowed to be used in 
a particular scheduling period, i.e., 

w{B{t)) = a ■ B{t), 0 < a < 1. (5) 

C. Energy Harvesting Model 

We assume a discretized model for the energy arrivals [9] where H{t) is modeled as an ergodic Bernoulli 
process (which is a particular case of a Markov chain). As a result, only two values of harvested energy are 
possible, i.e., H{t) G {0,e}, where e is the amount of Joules contained in an energy packet. The probability of 
receiving an energy harvesting packet during one scheduling period depends on the actual harvesting intensity 
(in the case of solar energy, it depends on the particular hour of the day) and is denoted by p{t). Note that a 
higher value of p{t) will be obtained in scheduling periods where the harvesting intensity is higher, e.g., during 
solar presence such as during the day, and a lower value of p{t) will be obtained during periods of solar absence, 
such as during the night. 

D. System Assumptions 

Let us collect all the channel gains, hk that includes the antenna gains, the path loss, and the fading, in h = 
{hk, yk G /GyU/Ci)}. Generally, the wireless channels depend on the specific scheduling period, h(f), as fhey vary 
over lime bul for simplicity in the notation, we will just refer to them as h throughout the paper. The traffic power, 
PBs{t) from (1), can be splil into power for voice and data connections as PBs{t) = J^keiCv Ffc(h), 
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where Pj(h) and pfc(h) are the instantaneous powers corresponding to the transmission toward the j-th and fc-th 
voice and data user, respectively. Let PRAuit) = ^Bs(i) + ^cpich be the overall radiated power hy the BS. 

The set of voice users request a fixed data rate and we assume that just one WCDMA code is assigned to 
them. This is translated into a minimum signal to interference and noise ratio (SINR) requirement as follows: 


Mvpk{h))hk _^ 

0{PRADit) - pk{b))hk + 0-2 “ ^ 


yk £ ICv, 


( 6 ) 


where My is the spreading factor for voice codes, 0 is the orthogonality factor among DL codes [6], and 
is the noise power. For simplicity in the notation and tractahility, we will consider the following approximation 
throughout the paper: ^ 

OipRADit) - Pkih))hk + a‘^ ^ OPRADit)hk + cr^. (7) 


On the other hand, the set of data users request a flexible service rate. The instantaneous throughput in the 
wireless access channel achieved during one particular scheduling period hy the k-th user, rfc(h), is upper hounded 
hy the maximum achievable rate that the access network is able to provide, which is formulated as 


W 

rkO^) < nfc(h)—log 2 { 1 + 


MDPkih)hk 


( 8 ) 


nk(h){9PRAuit)hk + 0-2 

where Md is the spreading factor for data codes, W is the chip rate, and nfe(h) is the number of codes assigned 
to user k. Notice that we have also approximated the denominator within the logarithm as in (7). 


III. Problem Formulation 

Let us introduce the following set of definitions: r = {rfc(h), VA: G ICd}, P — {Pfc(h), V/c G ICv}-, P — 
{pfc(h), yk G ICd}, n = {uk, (h) yk G ICd}- We formulate an optimization problem for the resource allocation 
strategy with backhaul and energy constraints to be executed at the beginning of each particular scheduling period, 
which involves finding the optimum resource allocation variables, r, p, p, and n that maximize the minimum of 
the expected throughputs (note that if a scheduling criterion different from the maximin approach is to be taken, 
problem (9) could be extended by just reformulating the objective function accordingly): 


maximize 
r.p, p,n, PradI*) 

subject to 


min Eh[rfc(h)] 
ksK- D 


Cl : 


Mvpk{h)hk 

OPRAD{t)hk + 0-2 


> r, yke ICv 


Vt € Vo 


^\ICd\ 

w ( 

C2, : rfc(h) < ^fc(h)— log2 f 1 + 


MDPkih)hk 


nk{h){0PRAD{t)hk + cj2) 


C4 : rJ ^ Pk{h) + pk{h) ] <<P{B{t)) 

\k£Kv 

C5 : nfc(h) < iV„ 


k&K.r 


k&K-L 


(76 : rfc(h) > 0, pk{h) > 0, nfe(h) >0, V/c G ICd 
C7 : PRAoit) = ^ Pfc(h) + ^ PA:(h) + Fcpich, 


k£lCv 


keKr 


yk G ICd 


( 9 ) 


'if the number of users is relatively high, PnAoit) ^ Pk, and the approximation is fair. In any case, the approximation provides a lower 


bound of the actual SINR value. 
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where > 1), is an overhead considered for the data transmissions to he sent through the hackhaul, Rbh{\KIv\) 
is the hackhaul capacity used hy the voice users^, being |/Cy| the number of voice users, Rbh is the overall 
backhaul capacity, F is the target SINR for the voice users, the function is related to g{-) in (4) as (f) = 

g{B{t)) — Tg ■ (-PcpiCH + Pc), and iVmax is the number of available codes for the data users. Although all the 
variables in the optimization problem (9) depend on the scheduling period t, we only keep such explicit dependence 
w.r.t. time in variable PRAD(f) to make explicit that that temporal evolution of the battery levels has a direct 
impact on the maximum power to be spent for the voice and data traffic, which is not constant along time. 

It is important to realize that problem (9) may not be feasible due to constraint Cl as it may happen that 
there could not be enough power to satisfy all the target SINR simultaneously. However, let us consider initially 
through the development that the problem is feasible (the feasibility condition will be developed later on the 
paper). Notice that, at the optimum, C4 is attained with equality. Otherwise, if C4 is not fulfilled wifh equalify, 
we could re-scale all fhe powers wifh a common positive facfor higher fhan 1 unfil C4 is fulfilled wifh equalify. 
This would increase fhe objecfive funcfion and all fhe ofher consfrainfs would sfill be fulfilled. Because of fhis 
reason, we can assume fhaf fhe optimum value of PRAD ( t ) is ^rad(^) = + Pcpich and we can eliminafe 

consfrainf C7 from problem (9). Consfrainf (72 sfafes fhaf fhe average fhroughpuf fhaf a user is experiencing in 
fhe access nefwork should nof exceed fhe maximum backhaul rale assigned lo such user (every user has been 
already assigned a porlion of fhe backhaul, as commented before). Again we slate fhaf Ibis is so because, in 
general, backhauls are nof able lo differenliale among users. If fhe backhaul was able lo do so, Ihen (72 could 
be rewriffen as . in any case, notice fhaf (he inslanlaneous rates allocated (o 

one user in (he access nefwork can be higher in some scheduling periods fhan (he maximum backhaul per-user 
rale thanks to the fact that queues are considered at the entrance of the access network. The 

average rate constraint (72 assures that the queues will be stable. 

It is easy to realize that the problem is separable into voice and data users without loss of optimality. For this 
reason, we start by analyzing the voice users. 


A. Resource Allocation for Voice Users 

Voice users must satisfy a minimum SINR constraint that is related to the target data rate service. Such 

Mvpk{^)hk 


constraints are: 


> r, VA: G JCv- 


( 10 ) 


dPlAu{'t)hk + 0-2 

It is straightforward to obtain the optimum power allocation for each voice user as follows (realizing that at the 
optimum, constraints (71 are fulfilled with equality): 

' Myilk 

At this point, we could check the feasibility of (9). The problem is feasible if 


( 11 ) 


Tg Y. Plih) < f{B{t)), 


( 12 ) 


keiCv 


which could also be written only in terms of the channels of the voice users, the current battery level, and some 
constants as follows: 


^ < Ki<p{B{t)) - K2, 




hk 


(13) 


^The overall backhaul capacity required for a set of voice users to be sent through the backhaul generally depends on the current number 
of voice users being served. In real deployments, voice users can be jointly encoded and, thus, the overall overhead for voice users may 
be reduced as the number of voice users increases. Anyway, in the problem formulation and the for the sake of generality, we just use 
the notation RBH{\iCv\)- 
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where ki = and K 2 = if the problem is not feasible, then we should consider either 

reducing such minimum SINR requirements (which will increase the constant term ki), or drop out some voice 
users in the scheduling period, or increase (j){B{t)) by increasing the value of a, but always guaranteeing that 
the maximum radiated constraint is not exceeded. 


B. Resource Allocation for Data Users 

Now that we have considered the voice users, we can tackle the resource allocation problem for the data users 
by solving problem (9). Note that problem (9) is convex once we know P^^it). To solve problem, (9), we will 
reformulate it by introducing the slack variable s, which preserves convexity [10], as 


maximize s (14) 

s, r, p, n 

subject to C2 ,..., (76 of problem (9) 


C8 : s < Eh[rfe(h)], \fk G )Cd 
Rb h - Rbh{\K^v\) 


(79 : 0 < s < 


Notice that we have introduced an additional constraint, (79. As it is clear from the formulation, this constraint 
does not affect the optimum solution, but it will help in the numerical search of the optimum value of the 
new slack variable s. Notice also that the previous optimization problem is time-coupled (we require the future 
channel realizations due to the expectation operator appearing in (78). In order to deal with such difficult problem 
involving expectations, we propose to use a stochastic approximation that has been proposed in the literature [11]. 
In this approach, the constraints involving expectations are dualized, and their Lagrange multipliers are estimated 
stochastically at each period. 

Let us start by dualizing constraint (78. Let A = {A^, \/k G ICd} be the vector of Lagrange multipliers 
associated to (78. The partial Lagrangian is given by Ccsis, A) = —s + “ ®"h[?’A:(h)]). In order to 

find the optimum s we have to perform the following minimization: 


minimize Ccs (s; A). 

C\^ 

-- CI/CdI 


(15) 


According to author in [11], when the objective function is linear in the optimization variable, the stochastic 
primal-dual algorithms present some numerical problems. This can be avoided by transforming the objective 
function introducing a general differentiable monotonically increasing cost function U{-) (e.g., the logarithm). 
Note that the introduction of such function does not modify the optimal value of the optimization variables (i.e., 
the solution is the same). Given that, setting the gradient to zero, Vs£c'8(S)'^) = 0 and solving yields: 




f (l^v I) 


s*(A) = {U) 


-1 






(16) 


\kCiKL 


where [/(•) is the derivative of [/(•) and (U)~^{-) is the inverse function of Uf). Once we know the optimum 
s*, the problem (14) is updated as follows (where we have skipped in the objective function the term that does 
not depend on the optimization variables remaining in the optimization problem): 


maximize 

r, p, n 


subject to 


Y. ^fcEh[r’fe(h)] 

(72,..., (76 of problem (9). 


( 17 ) 
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Now, we proceed to dualize constraint (72. Let /r = \/k G /Cd} be the vector of Lagrange multipliers 

associated to (72. The partial Lagrangian is 


^C2(r'fc(h), A,/i) = - ^ AfcEh[rfc(h)] + ^ //fc I Eh[rfc(h)] 


k&K 
= - Eh 


k&KL 


Rbh - Rbh{\K^v\)\ 


(Afc - /Ufc)rfc(h) 


.fceK:£ 


k'k 

ksK. D 


/ Rbh - Rbh{\K^v\)\ 


J’ 


V 






(18) 


(19) 


For given Lagrange multipliers A and /r, the optimization problem (14) is equivalently reformulated as (where 
we have skipped again in the ohjective function the term that does not depend on the optimization variables 
remaining in the optimization problem): 


maximize 

r, p, n 


subject to 


^ (Afc - Atfc)r’fc(h) 

(73,..., (76 of problem (9). 


( 20 ) 


Notice that the expectations are no longer present in the formulation because the remaining constraints (73—(76 are 
applied to instantaneous resource allocation variables (without expectations) and also because the maximization 
of the expected value of the objective function with respect to r, p, and n in the current scheduling period, in 
this case, is the same as the maximization of the term within the expectation. The problem now resides in the 
computation of the optimum Lagrange multipliers which requires knowing the statistics of rfc(h). If we solve 
the dual problem of (20), i.e., sup_;^^Q inf £(r, p, n, A,/i), where ^ means element-wise inequality and the 
Lagrangian, C, is defined in Appendix A, using a gradient approach, then the optimum multipliers could be found 
as [10]: 


A<»+">= (a<«> + e («*(AW) - El [rj (h; A<«>.m“)]))“ , Vt. 






Eh 




Rbh - Rbh{\ICv\)\ 




J 


yk, 


( 21 ) 

( 22 ) 


where e is the step size. Note that, it is not possible to compute the value of the Lagrange multipliers in real 
time, and then solve (20), as they depend on the statistics of rfc(h) that is a function not known a priori (it is the 
solution of the optimization problem itself). In this situation, we propose to follow a stochastic approximation [11] 
and eliminate such uncertainty constraint by estimating the multipliers stochastically at each scheduling period 
(with a noisy instantaneous unbiased estimate of the gradient) as follows (note that this philosophy is similar to 
the instantaneous estimation of the gradient in the LMS algorithm [12]): 


Xk{t + 1) = (Afc(f) + e (s*(A(f)) - rl{h; X{t), /r(f))))“ , V/c, 

Rbh - Rbh{\R'v\)\ 


Hk{t + 1) = Mfc(f) + € rl{h; X{t), n{t)) - 




J 


yk. 


(23) 

(24) 


The advantages of the stochastic techniques are threefold: i) the computational complexity of the stochastic 
technique is significantly lower than that of their off-line counterparts; ii) stochastic approaches can deal with 
non-stationarity environments; Hi) the distribution of the involved random variables h is not required. 

Once we update the values of the Lagrange multipliers, problem (20) can be solved using, for example a 
primal-dual approach. Notice that constraint (73 can be put directly in the objective function as, at the optimum. 


it is fulfilled with equality, i.e., r*(h; A(f),/r(f)) = n*(h; A(f),/r(f))^ loga [l + nUh;A(t),M(t))(eP,v(0fe.+^^) 

Thus, the resource allocation problem to be solved at the beginning of the scheduling period t is given by 




E w 

(Afc(f) - Hk{t))nkY;^ log2 

k^rC D 

subject to (74,..., (76 of problem (9). 


1 + 


MoPkhk 


nk{OPlpj^{t)hk + 


(25) 



Notice that we have not considered the dependency of the optimization variables with respect to the channels 
explicitly for the sake of simplicity in the notation. Problem (25) can be solved using, for example, a primal-dual 
approach, as it is given in Appendix A. 

It can be shown that the sample average of the stochastic rates, r^(h; A(f),/i(f)), satisfies all the constraints 
in (9) and incurs minimal performance loss relative to the optimal (off-line) solution of (9). This can be stated 
rigorously as follows: define F{t) = min^g^:^ I X]t=i f* minimum value of the 

objective function in (9). Then, it holds with probability one that as f — )■ oo: /) the solution is feasible; and ii) 
F{t) < f* -f ()(e), where (5(e) —)■ 0 as e —)■ 0. A proof of this result is not presented here due to space limitations 
but it can be derived following [11]. Let us introduce some important remarks here regarding problem (25): 

Remark 1: the values of {A^} measure how far the average rate s*(A(f)) is from the instantaneous rates 
served to the users in the access network. If the quality of the channels or the available powers are such that the 
instantaneous rates served in the access network are far from the target average rate s*(A(f)) for all the users, 
the sum of {A^} will increase and the system will reduce the target average rate s*{X{t)). 

Remark 2: in the access network, the values of {A^} are in charge of ensuring that all average rates tend to grow 
simultaneously (maximin objective function). Note that, if at any temporal period, s*(A(f)) 
then Xk{t + 1) grows and the priority of the k-th user to be served increases. 

Remark 3: if at any temporal period, r^(h; X{t), then fik{t -|- 1) > Hkit). For a fixed 

set of {Afc}, if Xk{t + 1) — fik{t + 1) decreases, the user will have a lower priority to be served in the next period. 
The same reasoning could be applied if r^(h; A(f),//(f)) < ^^h-Rbh{\iCv\) reverse situation. 

C. Resource Allocation Algorithm 

In this subsection, we just present the overall algorithm to solve the resource allocation for the voice and 
data users based on the approaches presented in previous sections. The stochastic updates, as well as the battery 
update, are also presented. The algorithm is summarized in Table I. Notice that, in the algorithm, steps 7 to 15 
correspond to the steps presented in Appendix A to solve the convex optimization problem (25). 

IV. Numerical Evaluation 

In this section we evaluate the performance of the proposed strategy. The scenario under consideration is 
composed of 1 BS, and 3 voice users and 6 data users. The maximum radiated power is = 9 dBm, the 

pilot power is Pcpich = 4 dBm (which represents the 13% of the maximum radiated power, as we considered 
in [13]), and the fixed power is Pc = 3 dBm (considering the model in [8], which was applied in [13]). The 
number of available codes for data transmission services is A^max = 15. All the users are mobile with a speed of 
3 m/s. The instantaneous channel gain, hk, incorporates antenna gains, Rayleigh fading with unitary power, and 
a real path loss of a town in Peru known as San Juan (see details in [13]). The orthogonality factor is 0 = 0.35. 
The code gain of data codes M/j = 16 and the minimum SINR normalized with code gain for voice users is, 
= —13.7 dB which corresponds to a rate of 12.2 Kbps. The noise power is a'^ = —102 dBm. The battery 
capacity is = 410 /iJ, the energy packet size is e = 30 fiJ, and a = 0.3 unless otherwise stated. The 

scheduling period for the data users and voice users are 2ms and 20ms, respectively, thus, Tg = 2ms. The utility 
function is U{-) = log(-). Two backhaul capacities have been considered in the simulations: Rbh = 2 Mbps and 
Rbh = 500 Kbps. The amount of backhaul capacity required by the 3 voice users considered in this deployment 
is RehUP-vI) = 173 Kbps. The overhead for the data transmissions is .^ = 1.2. The step size for the update of 
the stochastic multipliers is e = 10“^. For a more detailed description of the simulation parameters see [13]. 

In the simulations, we consider as a benchmark the case where the BS is connected to the electric grid (which 
means equivalently that the battery remains full of energy for the whole simulation). For comparison purposes, we 
also show the resource allocation of the proportional fair (PF) strategy [ 14] with an instantaneous per-user backhaul 



TABLE I 

Algorithm for Solving Resource Allocation Problem (9) 
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1: initialize Xk{t) > 0, > 0, Vfc € /Cd 

2: compute Prad(^) = + Pcpich 

3: Voice users 

4: compute pj^(h) = ^ yj. g 

5; if Ts ^ drop some voice users or reduce F, then go to 4 

6: Data users 

7: repeat 

8: initialize n ^ 0 

9: repeat 

10: =pI /?(«), A(t),/r(t)) using (27), Vfc G ICd 

11: A(t),/r(t)) using fixed-point iteration in (28), V/c G JCd 

12: until and converge 

13: update the dual variables, a^d using p^^^ and with (30) and (31) 

14: until and converge 

15: compute rj(h; A(t),/r(t)) withp^(n,/3) and nj(p,n, (p) 

16: update (dualized) primal variable: 

^bh-^bh(\^v\) 

17: s*{X{t)) = (((7)-i 

18: update stochastic dual variables: 

19: \k{t + 1) = (A(t) -P e (s*(A(t)) - r^(h; A(t), /x(t))))“ 

20: pk{t + 1) = (mW + e (r^(h; A(t), /x(t)) - J 

21: update battery with consumed energy and harvesting: 

Dmax 

22: B{t+l) = {B{t)-E{t) + H{t))^ 

23: t i — t -P 1 and go to 2 


constraint, rfc(h) < ^bh-R^h{\k.v\) ^ instantaneous sum constraint, J2k&K.n ■ Note 

that our proposed stochastic scheduling strategy resulting from (9) cannot he compared directly in a fair way 
with the PF approach under a sum-rate constraint. The reason is that in (9) we divide equally the hackhaul 
capacity among users in average terms whereas in the PF with sum-rate constraint, the overall hackhaul capacity 
is not forced to he distributed equally. The effective length of the exponential window in the PF has been set to 
Tc = 500 [ 14 ] 3 . 

Fig. 2 presents the instantaneous data rates served at the access network of four data users out of the six. In 
this case, the BS is connected to the electric power grid. As we can see, the instantaneous rates are able to exceed 
the backhaul capacity in particular scheduling periods whereas, at the same time, the average rates fulfills the 
maximum backhaul capacity as it is shown in Fig. 3. 

Fig. 3 shows the time evolution of the expected data rates of the three approaches. At any time instant the 
expected rates have been estimated using rk{t) = \ Et=i ^k{T)- We also plot the time evolution of s*{X{t)) and 
the per-data user backhaul rate. In this case, the BS is connected to the electric grid. The backhaul capacity is 
Rbh = 2 Mbps. Initially, we assume that the queues at the access network are sufficiently full so that all the bits 
demanded by the users are served. This makes the initial average rates violate the backhaul capacity constraint 

^The weights of the PF scheduler are calculated as uj^it) = where Tk{t) is the average throughput of user k computed as 

Tk{t) = (1 - ^jTk{t - 1) -f 
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for a short period of time (see the initial transient in the figure). This is due to the stochastic approximation of 
the multipliers hut, in any case, when the average rates converge, they fulfill all the constraints of the original 
problem. As we can also see from the figure, the limitation of the rates comes from the limited resources 
available at the access network, i.e., the power and the codes, as the backhaul capacity is not reached. It should 
be also emphasized that, the proposed stochastic approach provides a solution that introduces more fairness when 
compared with the PF approach as the average rates for the different users are quite similar. Fig. 4 depicts the 
same curves but now considering a backhaul capacity of Rbh = 500 Kbps. As we can see, now the system is 
limited by the backhaul and not by the limited resources of the access network. 

Fig. 5 shows the evolution of the stochastic estimation of the Lagrange multipliers, i.e., Xk{t) and /Ufc(f), for 
the cases where we have a backhaul capacity of Rbh = 2 Mbps and the case of having a backhaul capacity of 
Rbh = 500 Kbps. From duality theory, we know that if the backhaul constraint is not active, i.e., if the expected 
rates are below the backhaul capacity, then the optimum value of the multipliers is zero. This is what we see 
in the figure for the case of having a backhaul capacity of 2 Mbps. On the other hand, if the system is limited 
by the backhaul, then the optimum Lagrange multipliers are generally not zero as the corresponding constraints 
become active. From the figure we see that the multipliers converge to a non-zero value. The convergence of the 
multipliers states that the overall stochastic approach is working properly. 

Fig. 6 depicts the sum of the average data rates Ylk&Ko ^ function of the overall backhaul 

capacity {Rbh) for the different approaches when the BS is connected to the grid and when the BS has a finite 
battery with different harvesting intensities p. The black dashed line shows the total available backhaul for data 
users. In this figure, we are able to identify when the system is limited by backhaul and when the system is 
limited by the access network. For example, for the stochastic case connected to the grid, the system is limited 
by the access network when the capacity of the backhaul is above 1.2 Mbps. Concerning the comparison between 
our strategy and the PF approaches, we would like to emphasize again that the comparison with the PF approach 
under sum-rate constraint is not fair since in the last case the rate in the backhaul is not forced to be equally 
divided among users. Note, however, that even in that case, the sum-rate achieved by the PF is not much higher 
than the one achieved by our proposed stochastic approach. In addition, our proposed solution provides much 
more fairness as shown in Fig. 3 and Fig. 4. 

Finally, Fig. 7 depicts the time evolution of the instantaneous battery level of the BS when the electric grid is 
not available for the stochastic approach. We assume that the probability of receiving an energy packet during one 
scheduling period is p = 0.4 and p = 0.8 and two values of a have been considered: a = 0.1 and a = 1. Recall 
that a = 1 means that all battery could be used during one particular scheduling period (if the physical limitation 
of the BS allows it). In can be proved that, if no physical power limitation exists at the BS and the battery never 
reaches its maximum value, then, theoretically the expected value of the battery is given by 6 = 

However, because of the maximum power radiation at the BS and the battery overflows, the previous expression 
yields a lower bound of the true expected battery level, i.e., b < limt_ 5 ,oo IE[7?(f)]- For example, if p = 0.8, 
a = 0.1, then b = 240 pJ, but the figures shows limt_^oo IE[H(f)] = 340 pJ. The larger the value of a, the looser 
the lower bound is. 

V. Conclusions 

In this paper, we have proposed a resource allocation strategy based on the maximization of the minimum 
average data rates in a WCDMA system. By the use of stochastic optimization tools, we are able to consider a 
backhaul capacity constraint in terms of the average data rate, allowing the access network to offer higher rates 
taken advantage of good instantaneous wireless channel conditions. We have assumed that the BS is powered 
with a finite battery that is able to be recharged by means of a harvesting source. The dynamics of the energy 
harvesting, the energy spending, and the battery have also been taken into account explicitly in the proposed 
resource allocation problem. Simulations results showed that the proposed approach achieves more fairness among 
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the users when compared to the traditional PF strategy and, for some hackhaul capacities, the sum-rate is higher 
if compared with the PF with an instantaneous per-user hackhaul constraint, or the same if compared with the 
PF with an instantaneous sum-rate hackhaul constraint. 


Appendix A 


In this appendix, we present the technique to solve (25) assuming that the Lagrange multipliers are known 
(therefore, we omit the explicit dependence of the optimization variables with respect to the stochastic Lagrange 
multipliers, \{t), n{t), due to simplicity in the notation). Let /3 and (p he the Lagrange multipliers associated 
to constraints (74 and (75. There is no need to dualize constraint (76 because the solution will turn out to 
automatically satisfy it. The Lagrangian of problem (25) is given by 

MnPkhk 


^ (P, n,/ 3 , (/?) =- 


EW 

{Xk{t) - Pk{t))nk—log 2 

rvG/C D 


1 + 



nk{0P^f^^{t)hk + 0-2 

Ilfc -^max 


(26) 



For given Lagrange multipliers, /3 and p, we need to minimize the Lagrangian w.r.t. the primal variables. As it 
will be shown next, the structure of C (p, n,/), (/9) allows the minimization w.r.t. p and n to be found in closed- 
form. Because C (p, n, /3, p) is strictly convex and differentiable w.r.t. p and n, minimization w.r.t. these variables 
requires to equating the corresponding partial derivatives of C (p, n, /?, p) to zero. Differentiating the Lagrangian 
w.r.t. the data powers, equating the derivative to zero and solving such expression for the data powers yields 

/ {\k{t) - pk{t))nkW _ Uk {0PRADii)hk + 0 -^) 

I ln(2)/3MD Mnhk 

\ / u 

where the projection on the nonnegative orthant guarantees that constraint (76 is fulfilled. Proceeding similar with 
the optimum code allocation, we set the partial derivative of the Lagrangian w.r.t. Uk, equate it to zero and solve 
such equation for the codes, yielding: 

V-l 


Pfc(n,/3) = 


(27) 


n. 


(p, !!,(/?) = 


(Afc(f) - pk{t))WMoPkhk {{OP^pj){t)hk + 0-2) ln(2))’ 


MoPkhk 


GPRADit)hk + 


(28) 


(Xkit) - Pk{t))Wlog (l + 

where, also in this case, the projection on the nonnegative orthant guarantees that constraint (76 is fulfilled. Notice 
that a fixed-point iteration to compute the optimum code allocation, n\{p), can be used in this case. 

Having obtained the optimum primal variables as a function of the Lagrange multipliers, we now seek to 
find the optimum Lagrange multipliers to obtain the global optimum primal variables. The approach we propose 
to find the optimum multipliers is based on the maximization of the dual function D{/3,p) [10], which is 
defined as the minimization of the Lagrangian w.r.t. the primal variables, i.e., D{/3,p) = infp^„ £ (p, n, 
>C(K(n, /?)}, {n^(p, n, (/?)},/), (^). Then, the the multipliers are obtained by solving the dual problem as 


maximize D{/3,p) (29) 

subject to /3 > 0, (/^ > 0. 


Recall that the dual problem is always convex w.r.t. the dual variables and, thus, can be efficiently solved 
with a projected gradient method (if D{/3, p) is differentiable) or a projected supergradient method if it is not 
differentiable [15]. A valid supergradient for each particular dual variable is given by the constraint it is associated 
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with [15]. The update equations are given hy 

rf - {- E 11 . PO) 

V Kk&K-D V ® k&K-v / / / 0 

( Y1 4 "^ - ) , (31) 

V Vfce^CD / / 0 

where q indicates the iteration and the step size defined as (||VD|| 2 )~^, being VD the overall super¬ 

gradient of the dual function, is chosen such that the diminishing conditions are fulfilled, i.e., limfc_^oo = 0, 
Yhk&Ko ~ know the optimal dual variables, {3* and we can obtain the optimum power 

and code allocations, pXiP*) and The proposed iterative algorithm is based on the primal-dual block 

coordinate descent method for the update of the primal variables pk and [15]. 
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Fig. 1. Architecture of the target rural scenario under consideration in the paper. The BS is powered with a solar panel and a battery 
and the backhaul considered is based on WiFi-LD. The specific details of the real deployment as well as the location will be explained 
in the simulation section. 
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Fig. 2. Time evolution of the instantaneous data rates served at the access network and the backhaul capacity limitation per user with a 
backhaul capacity of 2 Mbps. 
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Fig. 3. Time evolution of the data rates for the different approaches and the hackhaul capacity per user when BS is connected to grid 
with a hackhaul capacity of 2 Mbps. 



Fig. 4. Time evolution of the data rates for the different approaches and the hackhaul capacity per user when BS is connected to grid 
with a backhaul capacity of 500 Kbps. 
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Fig. 5. Time evolution of the stochastic Lagrange multipliers for different backhaul capacities. 
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Fig. 6. Sum rates as a function of backhaul capacity for different approaches and different probability of energy packet p. 
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Fig. 7. Battery evolution of the proposed stochastic approach and the PF with sum constraint with a probability of energy packet p = 0.4 
and p = 0.8 a for a = 0.1 and a = 1. 

























































































































































































































